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Abstract: In this paper possible causes of various instability and chances of system break down in a 
power system network are investigated based on theory of nonlinear dynamics applied to a Power system 
network. Here a simple three bus power system model is used for the analysis. First the routes to chaotic 
oscillation through various oscillatory modes are completely determined. Then it is shown that chaotic 
oscillation eventually leads to system break-down characterized by collapse of system voltage and large 
deviation in Generator rotor angle (angle divergence), also known as chaos induced instability. It has 
been shown that chaos and chaos induced instability in Power system take place due to the variation in 
system parameters and the inherent nonlinear nature of the power system network. The relation between 
chaotic oscillation and various system instabilities are discussed here. Using the simple power system 
model, here it is shown that how chaos leads to voltage collapse and angle divergence, taken place 
simultaneously when the stability condition of the chaotic oscillation are broken. All nonlinear analysis 
is implemented using MAT LAB. It is indicated that there is a maximum lodability point after which the 
system enters into instability modes. All these studies are helpful to understand the mechanism of various 
instability modes and to find out effective anti -chaos strategies to prevent power system instability. 
Keywords: Angle divergence, bifurcation, chaos, Period Doubling, Power system, Voltage collapse. 



A power system is inherently of nonlinear nature i.e. the power system dynamics is described by a set 
of nonlinear equations obtained from system modeling and parameters. To a large extent, this is also due to the 
fact that most of the major power system breakdowns are caused by problems related to the system dynamic 
responses. It is believed that new types of instabilities emerge as the system approaches the stability limit 
which cannot be explain completely or accurately using linear analysis. So, theory of nonlinear dynamics is 
used to find out the proper explanation of those instabilities. Being an inherently nonlinear system, power 
system undergoes changes in state either qualitatively or quantitatively with the changes of one or more system 
parameters. This phenomenon is known as Bifurcation. Sometimes variation in parameter may result in 
complicated behavior which is random and indeterministic, known as chaos. This chaotic oscillation may lead 
to various instabilities. Most common types of system instabilities, which occur when the system is heavily 
loaded, are voltage collapse and angle divergence which eventually cause system breakdown. 

Voltage collapse in electric power systems has recently received significant attention by researchers. A 
number of physical mechanisms have been identified which possibly leading to voltage collapse. In several 
papers [4-10] voltage collapse was viewed as an instability which coincides with the disappearance of the 
steady state operating point as a system parameter, such as a reactive power demand is quasistatically varied, 
which is known as fold or saddle node bifurcation of the nominal equilibrium point. Dobson and Chiang [1] 
first studied and analyzed voltage collapse, and concluded that this phenomenon occurs at a hypothesized static 
bifurcation of equilibrium points taking place as system loading is increased. The static bifurcation mechanism 
for voltage collapse postulated in [1] was investigated in [2] and in [3]. It is therefore not surprising that saddle 
node bifurcation is being studied as a possible route to voltage collapse [4]. 

Another possibility is that steady state operating point loses stability before the saddle node 
bifurcation. If this occurs in a given system, stability of the nominal equilibrium point may be lost prior to 
static/saddle node bifurcation point through a Hopf bifurcation. Study of Hopf bifurcation in power system was 
done in details in [2, 3, 4, and 5]. Except these, other type's bifurcation occurs in power system like Torus 
bifurcation [4], cyclic fold bifurcation, period doubling bifurcation [3, 2]. 
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In this paper a complete and detail bifurcation analysis has been done which shows different 
behavioral changes (Bifurcation) with the slow and gradual variation of load reactive power. It has been shown 
here how the stable oscillatory behavior of the power system model tends to chaotic instability through period 
doubling bifurcation. Period doubling bifurcation (PDB) is the most important route to chaos in power system 
model which is analyzed here with great emphasis. Except continuation method, a detail and explicit picture of 
PDB has been developed which clearly shows the internal behavior changes of the proposed system which 
eventually leads to chaos. This paper, for the first time, proposes that Voltage collapse and angle divergence 
phenomenon which make the system unstable, occurs simultaneously. Also this paper gives an indication on 
maximum lodability point after which system tends towards instability. 

The numerical solution of the basic nonlinear differential equations of the proposed model of power 
system network and load are implemented using the MATLAB environment [16] and assembly language 
programming. 

II. Brief Review Of Nonlinear Theory [17,18] 

A typical nonlinear system with state x can often be expressed as- 

x = f(x,fi);x E9t n , fi E9t p (2.1) 

The corresponding properties of such a system are: 

• The solution of (2.1) is called trajectory. With initial condition x(t 0 ) = x 0 , the solution is given by x(t ) = 

A t (Xo). 

• Four steady state behaviors are associated with the nonlinear system. Equilibrium points, Periodic 
solutions, quasi periodic solution and chaos. 

• The equilibrium points mean the solutions of the nonlinear equation / (x, A)= 0. Alternatively an 
equilibrium point is a degenerate trajectory which stays in the equilibrium point for all time. It is 
asymptotically stable if all the Eigen values of its corresponding Jacobian matrix have negative real part. 
Power system is generally operate on a stable equilibrium point. A t (x*) is a periodic solution if A t (x*) = 
A t+T (x*) for all t and some minimal period T>0. It represents a limit cycle which is a self sustained and 
bounded oscillation, and, is stable or unstable depending upon its characteristics multiplier.. 

• Finally, chaos is a random, indeterministic phenomenon exhibits stable, bounded but aperiodic behavior. 
While equilibrium points are zero dimensional and periodic solutions are one-dimensional, chaos is more 
complex and having fractional dimension. 

2.1. Bifurcation Theory: 

Bifurcation Theory is used to interpret the way in which qualitative changes occur in the system as 
one or more parameters are varied. A power system is modeled in the typical form of a nonlinear dynamic 
system with state x: 

x = / ( x , u); xG« n , G5R P (2.1.1) 

|i represents the vector of the system parameters that can be varied during the analysis. 
Bifurcations in dynamical system mean "qualitative" changes of the asymptotic behavior of the system 
trajectory (2.1.1) which is obtained by varying the |i components. 

At a value of □ = Q the vector field / loses its structural stability, is called the Bifurcation Point and Q. ,the 
Bifurcation value. This simply means that the Phase portraits for □ < Q. and for □ > are different. In power 
system's nonlinear model there are following types of bifurcations taken place depending upon the Jacobian 
(3) of the system (2.1.1) - 

• Saddle-node bifurcation (SNB):- A saddle-node bifurcation is a local bifurcation in which two fixed 
point (or equilibrium) of a dynamical system collide and annihilate each other, at this point the Jacobian has a 
zero Eigen value and no other Eigen value with zero real part. 

• Hopf bifurcation (HB):- A Hopf-bifurcation is a local bifurcation in which a fixed point of a dynamic 
system loses stability as a pair of the complex conjugate Eigen values of the linearized system around the fixed 
point cross the imaginary axis of the complex plane. If has 3 a pair of complex conjugate Eigen values on the 

imaginary axis where all other Eigen values are off the imaginary axis, then Hopf Bifurcation results the 
emergence of a family of periodic solution in the vicinity of |i c . If the periodic solution is unstable then it is 
called "subcritical" and "supercritical" when stable. 

• Period Doubling Bifurcation (PDB):- During this kind bifurcation with the change in multiplier a 
new periodic solution or orbit emerges from the previous/existing solution with periodicity approx. Twice that 
of the previous one. If there is a sequence of such bifurcation which accumulate at a critical value \i = □ c then 
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period almost becomes infinite, which means we get a aperiodic but bounded solution to the system which is 
called "Chaos". Actually this is one special type of Hopf Bifurcation. 

III. Power System Model For Nonlinear Analysis 

For applying theory of nonlinear dynamics, a simple power system model [1] is considered here. This 
model is widely used for nonlinear behavior study of power system [2, 3, 4, and 14]. In this power system 
model generator is represented by classical model. Here the system is represented by a set of four ordinary non 
linear differential equation. In this paper, classical model is implemented for a three BUS power system 
network as follows. 

A simple 3 -BUS power system shown in fig.l and its equivalent circuit in fig.2, this model consists of 
an infinite bus on the left, a load bus on the center and a generator bus on the right. Y 0 Z(Dr))nij L ) and Y m 

Z(D ^ ) are the admittances of the transmission lines. One of the generator buses treated as slack bus and 

iii 

the other is described by swing equation. The concept of an infinite/slack bus refers to a particular node of the 
system with enough capacity to absorb any mismatch in the power balance equations. Thus, it can be 
considered as fictitious generator with constant voltage magnitude E 0 and phase EJ> (usually E 0 =1 and EJ> = 0). 
On the other hand, the generator has constant voltage magnitude E m but the angle varies according to the 
so-called swing equation. 




Figure 1: 3-bus power system model 



Y0/(-90-n2) 



V/d 



Ym((-em-ir2) 



e 



EO ZiO " C 




P-Q 



9 



Em/ dm 



Figure 2: equivalent circuit 



s ni =□ 

M tti = -d m D +E m Y m sin (QJ + E m V Y m sin (□□£}+■ _„)] 

r2 



So, the set of four dynamic equations is developed for this model as - 

(3.1) 
(3.2) 

k qw J= □ l^V □ V v □ Q> □ Q + Q (On, V) (3.3) 

Tkq W kq V 2 t' — kp W kq V 2 V z + (k pw k qv □ Vk pv )V + k qw [P (CU V) □ Pi] □ V fQ (Qn, □, V) DQo □ Q 
(3.4) 

Where, Qn = Generator rotor angle, □ = Generator load angle, □= Angular frequency, V= Magnitude of 
Generator load voltage. 

The load bus, with voltage magnitude V and phase Q consists of an induction motor, a generic load P- 
Q and a capacitor C. The dynamics of this part is derived from a power balance at the bus. Considering an 
empirical model for the induction motor [9] and a static load P-Q, the active and reactive power supplied to the 
load is - 
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P (On, □, V) =P 0 + k pw dD+ k pv (V + TdV) + Px (3.5) 
Q (Qn, □, V) = Q)+ k qw dD+ k qv V + k qv2 V 2 + Qx (3.6) 

Where T, k pw , k pv , k qw , k qv and k qv2 are constants of the motor, P 0 , Qo and Pi, Qi are the static active and 

reactive power drained by the motor and by the load P-Q, respectively. 

Here in this power system network the reactive power demand at the load bus Qi is chosen as 

bifurcation parameter. Therefore the model has the form- 

* = /(x, □) 

where x= f CU □ , D,V] T is the state vector and □ = [Q] T is the bifurcation parameter vector so that increasing 
Qi corresponds to increase in load reactive power. 

IV. Numerical Solution- Simulation & Result 

Here the numerical solution of the nonlinear Eqn. (3.1 - 3.4) has been carried out using MATLAB 
simulation, the result of which can be sub divided into following two parts - 
1) Bifurcation analysis and Chaos, 2) Chaos induced instability and system collapse. 

4.1. Bifurcation and chaos: 

In this paper continuous Bifurcation diagram is plotted using MATLAB continuation software 
'MATCONT' to indicate different types of Bifurcation mentioned in II along with the respective bifurcation 
values in the proposed Power system model described by Eqns. (3.1-3.6). 

Fig.3 shows the complete bifurcation diagram using the continuation method depicting the variation 
of system voltage V (p.u) with Qi. As the power system network is inherently nonlinear system it's dynamic 
equations are also nonlinear so it's solution exhibits different dynamic behavior with changes in system 
parameters, which is completely identified in the bifurcation diagram in Fig 3. 
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Figure 3: bifurcation diagram using continuous method showing various type of dynamic behavior 



Suppose the power system model describe in section 3.1 is operating at a stable equilibrium point. 
Now Qi is slowly increased while other parameter remains fixed. At each parameter step, the system eigen 
values are calculated. In the course of computing stationary branch, three critical points S1,S2 & S3 are 
detected as shown in the Fig.3, at which system changes it's stability. In Fig.3 the solid and dotted lines 
represent the stable and the unstable stationary trajectory. The stable trajectory becomes unstable at a 
subcritical Hopf bifurcation (SI) for Qi=10.868 and regain stability at supercritical Hopf bifurcation (S2) 
at Q i=l 1.407. The Eigen values of two critical points SI and S2 are calculated numerically as (0.0000 ± 
j3.4234, -125.2174, -18.2564) and (0.0000 ± j2.8957, -92.1458, -2.3652) respectively i.e. at SI and S2, a pair 
of complex conjugate Eigen values cross the imaginary axes. So the real part of the complex conjugate Eigen 
values of the system becomes positive through SI and negative through S2. After a short stable region, one real 
Eigen value becomes positive at S3. This last critical point S3 is called Saddle node bifurcation point where 
one real Eigen value becomes zero. After that system becomes completely unstable. 

At point SI and S2, the possibility of oscillatory dynamic behavior is investigated by applying Hopf 
bifurcation theory. At this two points periodic branches are calculated. The periodic branch emanating from SI 
up to Pdl is unstable because it is subcritical. Then at Pdl periodic limit cycle emanates and gains stability. 
With further increase of Qi period doubling bifurcation occurs at Pdl where the previous periodic solution is 
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bifurcate to a new periodic solution or orbit with periodicity approx. twice that of the previous one. If one keep 
tracing an old periodic orbits, it again passes through the unstable periodic region between Pdl and S2 and 
reaches S2, the second super critical Hopf bifurcation point. Now if we decrease the value of Ql, it reaches 
second period doubling point Pd2. the periodic orbits gain stability again in the supercritical region between 
Pd2 and S2. After Pd2 if we decrease the value of Qi, system enters into a series of period doubling 
bifurcation. All these limit cycles represent periodic oscillation in system behavior with different frequency. 

After a short stable region beyond S2 system become unstable at S3. If the reactive power of the load 
Qi is increased beyond this value, the system looses it's stability and leading to system collapse. 



4.2. Chaos via sequence of Period doubling: 

According to bifurcation theory one way to chaotic motion is through a sequence of period doubling 
bifurcations. Therefore, we give special attention to two Period doubling points Pdl and Pd2, in this case and 
tried to show explicitly what happened after these two points with the variation of Qi with the detail 
bifurcation diagram plotted in the neighborhood of these two points. First, consider the period doubling 
bifurcation point Pdl at the left side of main bifurcation diagram of Fig. 3. With the further increase of Qi a 
new periodic orbit with periodicity two is emerged, and then the periodic orbit of period four and so on. After a 
certain value of Qi the system finally reaches to a state which is completely indeterministic and having infinite 
or fractional no of period or become aperiodic but exhibits a bounded oscillation. Dynamic behavior of the 
power system model at this region is very complex and unpredictable and there is no obvious relation between 
cause and effect. At this region the power system network shows completely erratic and random behavior. This 
bounded but random oscillation is known as Chaos. It can be termed as Left side chaos (LSC) as it occurs left 
part of continuous Bifurcation diagram of Fig.3. This Left Side Period doubling Bifurcation leading to 
Chaos is shown explicitly in Fig.4. 
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Figure 4: bifurcation diagram showing left side period doubling bifurcation leading to chaos 



In the right hand side of the periodic branch, the same phenomenon takes place but here the state 
enters into period doubling bifurcation at Pd2 towards the chaos with the decrease of values of Qi from its 
value at Pd2. This chaos can be termed as Right side chaos (RSC). This right Side Period doubling 
Bifurcation leading to Chaos is shown explicitly in Fig.5. 

Bifurcation Diagram showing Period Doubling 
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Figure 5: bifurcation diagram showing right side period doubling bifurcation leading to chaos 
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Table I: values Qj (bifurcation parameter) at different bifurcation point 



Value of Ql 


Nature of Bifurcation 


10.946 


Subcritical Hopf Bifurcation (SI) 


11.407 


Supercritical Hopf Bifurcation (S2) 


10.870 


Period- 1 LHS 


10.880 


Period-2 LHS 


10.883 


Period-4 LHS 


11.392 


Period 1RHS 


11.388 


Period 2 RHS 


11.384 


Period 4 RHS 



4.3. Phase Plots and Time plots: 

Below the various phase plots are shown which describes the various states of the system dynamics 
with different values of Qi, during PDB leading to chaos through both Left side PDB and Right side PDB. 
Also the time plots are shown for different values of Qi. These plots are obtained from the numerical solution 
of the Eqn. (3.1.1-3.1.4) which are numerically integrated with initial condition (0.3, 1.5, 0.2, 0.97) and 
(0.315, 0.150, 0.150, 0.98) using Runge-Kutta method. Here Chaos is observed for Qi=10.894, called Left 
side chaos and Qi=11.383, called Right side chaos. 
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Figure 6: phase plots and time plots showing left side PDB leading to chaos 
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Figure 7: phase plots and time plots showing right side PDB leading to chaos 



4.5. Chaos induced instability and system collapse: 
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Though chaos exhibits aperiodic and random oscillatory behavior and completely indeterministic, it is 
a bounded and marginally stable region. But Chaos is very sensitive to initial condition and system parameter 
variation. Any small change to them can break their stable oscillation. Here we discuss what happen after 
stable chaotic oscillation is broken in power systems. It will be shown that chaos can lead to Voltage collapse 
and angle divergence simultaneously when value of Qi is increased beyond the value corresponds to Left side 
chaos and decreased below the value corresponds to Right side chaos, which makes a stable system into 
complete unstable and the system breakdown takes place. 

• Voltage collapse and Angle Divergence: -Many studies have observed this phenomenon [1-11] taking 
place individually. But in this paper, voltage collapse occurring along with Angle divergence simultaneously is 
reported for the first time in Power system dynamic stability study. During Voltage collapse the system voltage 
sharply decline to a very low value and possibly brings the blackouts and Angle Divergence is the 
phenomenon when the generator loosing synchronism i.e. rotor angle difference is much more than 2n. Here 
critical point at which voltage collapse and angle divergence phenomena take place simultaneously is 
Qi=10.890 after Left side chaos and Qi= 11.377 after Right side chaos, using same model, parameter values 
and same initial condition. These are shown below in Fig 8 & 9 where time plot of load Voltage and Generator 
angle are given for the above mentioned values of Qi.For both the values of Qi we get almost same diagram. 
From the Fig 8 & 9 it is seen that voltage collapse and angle divergence appears after chaos is broken. 




Figure 8: Voltage collapse at Qi=10.890 and Qi=11.377 Figure 9: angle divergence at Qi=10.890 and 

Q 1= 11.377 



V. Conclusion 

In this paper various nonlinear dynamical behavior of proposed power system model has been deeply 
studied using both "Matlab" and Assembly language Programming. Here cascaded period doubling bifurcation 
which is one of the most important routes to chaos and system instability has been observed in details. All 
previous studies used continuation method and compact continuation software's such as AUTO [15] for 
Bifurcation studies in power system network which predicted the occurrence of Period Doubling Bifurcation 
leading to chaotic phenomenon but were not able to show explicitly and elaborately, the structure of period 
doubling bifurcation leading to chaos. This limitation is overcome in this paper where bifurcation analysis is 
done using both MATLAB based continuation software MATCONT to get a continuation Bifurcation diagram 
and also get a expanded and elaborate diagram of Period doubling bifurcation which clearly shows various 
state changes of the proposed system with the variation of Qi which eventually leads to chaos. This is also 
described by various phase plots and time plots for different values of Qi.The relationship between chaos and 
major instability modes in power system, such as Voltage collapse and Angle Divergence, has thoroughly been 
observed. Here it is shown for the first time that Chaos can induce voltage collapse and angle divergence 
simultaneously when the stability condition of the chaotic oscillation are broken. So it can be conclude that 
chaos is an inter mediate stage of the instability incident when changes in system parameter causes system 
breakdown. As chaos is very sensitive to initial condition and system parameters, any variation of them can 
make chaos to be annihilated and breaks in into instability. In a real power system all system parameters are 
fluctuating with changing operating condition and disturbances. 

So chaos possibly exists in power system network as in prior stage of instability. When disturbance 
happens, power system comes into a transient stage. If the disturbance is small, HB may happen and stable 
oscillatory behavior follows. If the disturbance is prolonged, system may come into chaos. And, when the 
disturbance becomes larger, the chaos may be broken. Voltage collapse, angle instability or voltage collapse 
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and angle divergence simultaneously may happen. If the disturbance is very large, system may directly come 
into the above three instability conditions over the stages of HB, chaos and chaos breaking. Here it is indicated 
that when there is a large disturbances in system parameter leading to system collapse, chaos is very likely to 
be an intermediate transient stage between a stable and unstable region a power system. All this studies are 
helpful in understanding, how various instabilities due to the interaction between system parameter itself take 
place and the possible routes through which the system moving towards breakdown through various stable and 
unstable oscillatory modes It also indicated that there is a maximum lodability point after which the system 
enters into instability modes. All these studies are helpful to understand the mechanism of various instability 
modes and to find out effective find appropriate measure to prevent nonlinearity induced instability in power 
system. 
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